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Abstract 



The hadronic mass distribution in semileptonic S-meson decays can 
be used for extracting the charmless part and thus determining the 
l^ub/^cfel ratio. We take into account first-order perturbative as well 
as non-perturbative QCD corrections. The sensitivity to model as- 
sumptions is studied and an estimate of the remaining uncertainties 
is performed. 



1 Introduction 



The CKM matrix element V u b plays an important role in the determination 
of the unitarity triangle. The cleanest method to obtain the absolute value 
\V u b\ is through the measurement of semileptonic b — > u transitions, which 
will eventually give the length of one of the sides of the unitarity triangle. 

The possibilities of determining V u b from semileptonic decays have been stud- 
ied in detail in the BaBar Workshop [1]. Here the conclusion was reached 
that the determination of V u b will be performed by a mixture of different 
methods. 

One of the theoretically cleanest possibilities is to use inclusive semileptonic 
decays. Placing a cut on the hadronic invariant mass of the final state can in 
principle eliminate the charm contribution, which otherwise would be over- 
whelming. This cut on the hadronic invariant mass can be implemented at 
the asymmetric B factories, making this method experimentally feasible. 

From the theoretical side the decay rate, including cuts on the lepton energy 
as well as on the hadronic invariant mass of the final state, can be computed 
systematically within the framework of the l/m& expansion, except for certain 
regions of phase space, where the l/m& expansion has to be replaced by an 
expansion in twist. To describe these regions of phases space one has to 
introduce a so-called "shape function" , which in principle introduces a large 
hadronic uncertainty. 

Another method that has been proposed [2] avoids the twist expansion and 
relies only on a standard 1/toq expansion. This method needs a measurement 
of the lepton invariant mass spectrum, in which the regions of phase space 
where the shape function plays a role are kinematically suppressed. This 
method will also allow a clean determination of V u b- 

An alternative approach has been put forth in [3, 4] where the factorization 
into soft, jet and hard sub-processes [5] has been employed to relate the 
radiative b decays to the semileptonic ones in a way which explicitly reduces 
the impact of the shape function uncertainties. Then a prediction of the ratio 
iKifc/VLl can be obtained with a good accuracy in a model independent way. 
This approach is similar to the one in [6]. 

Exclusive decays will open a completely different window on V u b, however, 
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in these decays a certain model dependence seems to be unavoidable, unless 
lattice data become reasonably precise. 

From this variety of methods to determine |14&|/|Kb|, we shall expand in this 
paper on the one in which a cut is applied on the hadronic invariant mass 
in semileptonic B decays to filter out the b — ^ u transitions. The advantage 
is that the hadronic invariant mass may be easier to measure, however, this 
method involves the shape function and potentially has larger theoretical 
uncertainties than the inclusive method using the leptonic invariant mass. The 
method based on the hadronic invariant mass spectrum has already been 
discussed in [7], where the main focus was on the perturbative contributions 
of order a s and flo a t- 

The purpose of this paper is to consider this approach in detail and to try to 
estimate the uncertainties, including perturbative as well as non-perturbative 
contributions and, in addition, a cut on the lepton energy. It turns out 
that the main uncertainties originate from the heavy quark mass m& (or, 
equivalently, from A = Mb — where Mb is the 5-meson mass) and the 
strong coupling a s , while the uncertainties introduced by the shape function 
are small. 

The next section deals with the kinematics. Then, in Section 3, we give the 
radiative corrections for the partonic process b — > uti>t to order a s . In sec- 
tion 4 we include the leading-twist non-perturbative effects; this requires the 
introduction of the light-cone distribution function for the heavy quark, for 
which we use a simple parametrization. We combine the perturbative and 
non-perturbative corrections and study the uncertainties in the determina- 
tion of |V^&|/|Kfi| in section 5. 

2 Kinematical Relations and Definitions 

We shall first define the kinematic variables for the partonic process b — > uli>i. 
Although we also use results for the semileptonic b —> c decay rate, the latter 
is considered in the partonic framework solely. In fact, we only need the total 
rate with a single lower cut on the electron energy for b — > c; we thus refer 
the reader to [8] for further discussion of the kinematics of this process. The 
initial state b quark has a momentum = m^v, where v is the velocity of 
the B meson. With the momentum transfer to the leptons q = k + k' the 
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variable p = q — m b v is the partonic momentum of the final state. Writing 
Eg = vk for the energy of the lepton, E p = vp for the partonic energy of 
the final state and p 2 for the partonic invariant mass, we define the following 
re-scaled variables 

_2E i _2E IL _l_ _l_ 

x , x p , y „ , z 2 . [i) 

m b m b ml m% 

One of these variables is superfluous when calculating the triple differential 
decay rate and can be substituted by using the relation 

x p + y - z = 1 . (2) 

Using these definitions, we compute the triple differential partonic rate to 
order a s and order l/m b , and split it into the tree- level term T° and the 
0(a s ) correction T 1 : 

37T 

The radiative corrections to 0(a s ) have recently been calculated by [9]; we 
have checked their result and find full agreement with ours. The formulae 
up to terms 0(a s ) are quite tedious and can be found in [9]. However, we 
present them in the Appendix in a form suitable for our numerical evaluation. 

We shall now discuss the hadronic kinematics. The momentum of the initial 
B meson is Mbv = (mb + A)v, where we have used the relation between the 
B meson and the 6-quark mass 

M B = m b + A (4) 

which holds to leading order in the l/m& expansion. Consequently, the 
hadronic mass of the final state is M\ = (Mbv — q) 2 = (p + Av ) 2 = 
p 2 + 2E P A + A 2 and hence involves both the partonic invariant mass and 
the partonic energy. Thus we have 

^YParton l f 2 ~ x Zm f ax j3pp er t 



dMldEi 01 J J p J dxdx p dz 

X 1 l-x " P 



(m&) = dx dx p / dz— — — — x 



5(M 2 X - m\z - x p Am b - A 2 )5 ^ - ^) . (5) 
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where the kinematic limits of the z integration depend on the cuts as well as 
on x p . They are 

_ _ / for x p < 1 , 

Zmin -\x p -l for x p > 1 , {0) 

and 

Z m ax = (1 - X)(x p + X - 1) . (7) 



3 Perturbative corrections 

The analysis we perform requires the knowledge of the triple differential 
partonic rate to order a s , as can already be seen from Eq. (5) as well as from 
the final formula, Eq. (12). The Born approximation is proportional to the 
delta function of argument z, 

= 12T (2-x-x p )(x + x p -l)5(z), (8) 



dx dx p dz 



where 



_ G 2 F m\ 

r ° " iW (9) 

As is well known, the virtual one-loop correction contains an infrared diver- 
gence that cancels with the real gluon emission. We shall discuss in the next 
section how to include the leading-twist non-perturbative effects by "smear- 
ing" over some small window in the hadronic invariant mass. While choosing 
an excessively small value for this window would still yield a meaningless 
result, a region of size Amj is believed to exist, where this smearing provides 
a realistic approximation. Instead of keeping track of the virtual and real 
parts, one can now use an appropriately integrated distribution, which is 
subjected to the smearing procedure. The tree-level term is rather simple to 
implement numerically but, on the other hand, we have found it very useful 
to perform one integration of the one-loop correction analytically. On in- 
spection of Eq. (12) in conjunction with Eq. (5), it is clear that the following 
integral is of use: 

T dxdxpdz' 
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This function is a lengthy and tedious expression, and it is given in the 
Appendix. 

4 Leading-twist non-perturbative corrections 

It has been shown that the leading-twist non-perturbative corrections can be 
implemented at tree level by redefining the heavy-quark mass and a subse- 
quent convolution with a so-called shape function [10, 11, 12, 13, 14]. This 
convolution corresponds to an integration over the light-cone variable k + , 
namely 

A 

dT hadron = J dT^™ {m b + k + ) f {k + ) dk + (11) 

in dT H ® ET . Although this formula is quite suggestive [15], it has been shown 
recently that it contains, in fact, spurious contributions of sub-leading twist 
[16]. Furthermore, once radiative corrections are implemented, this simple 
convolution formula will probably no longer hold [17]. However, there is no 
fully consistent way yet to include radiative corrections into this convolution, 
and hence we proceed in a naive way as suggested in [15]. 

This naive way is to actually use the convolution formula (11) also beyond 
tree level, which is at least as consistent as using the ACCMM model [18] 
beyond tree level, which is common practice. In fact, the connection between 
the ACCMM model and the shape function formalism has been pointed out 
in [19]. 

It is convenient to change the variable of integration according to 

M B 

dT hadron = J dT parton (m*) f (m*) dm* (12) 



with m* = nib + k + . 

The shape function is a nonperturbative function which has to be determined 
either from experiment or by some model. A few relations are known for the 
moments of the shape function 

(*°) = 1, (fci> = and (k 2 + ) = - 1 -X 1 (13) 
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For our study we shall use an ansatz for the shape function. Taking a simple 
three parameter function [20] 

f{x) = Ne cx {l-x)\ x = l- MB ~ m \ (14) 

we can express the parameters of this ansatz in terms of the HQET param- 
eters (13) 

c c 3A 2 
N = - — T , a = c-l, c= — — . (15) 
Ae c r(c)' Ai K J 

We shall discuss the dependence of our results on this ansatz in the next 
section. 



5 The measurement of \V u b\ 

For the measurement of |V^&| we propose to study semileptonic B decays with 
certain cuts. The first cut is on the lepton energy Ei, which is mainly given 
by the experimental limitations on detecting electrons with small momenta. 
The second cut is on the hadronic invariant mass M x of the final state, which 
serves to suppress charm. We define the semileptonic rates including cuts as 







2 M B 


J dM 2 x 







E C ut 



r(Ml t ,E cut )= I ml / dfl g^ggp* , (16) 



Clearly the region < Mx < is dominated by b — > u transitions, and we 
thus use in this region the expressions for b — > u decays only. 

We shall normalize everything to the rate with no cut on the hadronic in- 
variant mass, and thus obtain, for the ratio r: 



\V ub \\ T(Ml t ,E 



|y c6 | 2 r(Mg, E cut ) 

where in the denominator we can safely take into account the b — > c channel 
only, the b — > u contribution being only about 1%. 

For the charmless decay rate entering the numerator, we include the results 
discussed above. The perturbative corrections are taken into account to one 
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loop [8], while the non-perturbative ones are included to leading twist. The 
b — > c decay rate in the denominator is evaluated including 0(a s ) corrections. 

Unlike the approach advertised in [2] and [21], the shape of the hadronic 
invariant mass spectrum depends on the shape function, for which we use 
the parametrization (14). However, as it will turn out, the precise form 
of the shape function is irrelevant for the ratio (17), the main sources of 
uncertainties are the value of the strong coupling and that of the quark 
masses. The latter may be replaced by A, we thus use m h = Mb — A and 
m c = M D — A, where M denotes the spin- averaged meson mass. 

The dependence on the shape function can be studied by checking the sensi- 
tivity of our result with respect to higher moments of this function. Since we 
have not included any 1/m 2 , terms in our calculation the ratio (17) should be 
reasonably insensitive already to the second moment of the shape function, 
which is given in terms of the kinetic energy operator Ai in (13). This forces 
us to make M cut as large as possible. 

We have examined the ratio defined in Eq. (17), allowing the variations 



0.4 GeV < A < 0.75 GeV, 0.2 < a s < 0.3 and -0.6 GeV 2 < Ai < -0.1 GeV 2 . 



The main uncertainty is induced by a s and A, where A is equivalent to 
the heavy quark mass. It has been argued that these two quantities are 
correlated; the size of the radiative corrections depends on the particular 
choice of the mass. 

Using the pole mass scheme, it has been shown that the radiative corrections 
are large and the perturbation series converges very slowly. Treating both the 
pole mass and the perturbative contributions independently, the uncertain- 
ties of the quantities would simply add, leaving us with a large (and certainly 
overestimated) uncertainty. In this case we get (see the lighter band in Fig. 



where we have used ml° e = (4.75 ± 0.15) GeV and a s between 0.2 and 0.3. 
In this way we obtain a theoretical uncertainty of 25% in the determination 
of the ratio. 

Another option is to switch to a short-distance mass definition, such as m^f s 
by replacing 



1): 



1.15 < r < 1.96 at M 2 cut = 4 GeV 2 



(18) 




(19) 
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Figure 1: Ratio function of the hadronic mass cut, with (dark) or 

without (light) correlation between the strong coupling constant and the 
pole mass 



which reduces the size of the coefficients of the perturbation series, and the 
pertrubative uncertainties thus become smaller. Using a recent value for 
M^ s [21] 

mf s = (4.25 ± 0.08)GeV, (20) 

we arrive at an estimate for r with a smaller uncertainty (the darker band in 
Fig. 1) 

(21) 



1.28 < r < 1.85 at M 2 cut = 4 GeV 2 . 



In order to display the dependence on the input parameters, we choose the 
"average" values of the three parameters: 

Javer = q 55 ^aver = q 3^ yaver = _q 35 ^ 

and obtain up to linear terms in the variations 



AA Ao s AAi 
r = 1.58 - 0.86^— - 0.26 ^ + 0.06- 



A° 



a" 



Af 



(23) 



8 



This explicitly shows that the dependence on Ai, which is the second moment 
of the shape function, is weak. The dependence on even higher moments is 
expected to be further suppressed by inverse powers of the heavy quark mass. 

In conclusion, our best estimate for r is 



which corresponds to a theoretical uncertainty in the determination of \ V u b/V c b 
of about ten percent 



We have performed a detailed analysis of one of the possibilities to obtain 
V u b from inclusive semileptonic B decays by placing a cut on the hadronic 
invariant mass to get rid of the charm background. This method has been 
criticized since it depends on the shape function, which describes the endpoint 
of the hadronic invariant mass spectrum. This function is not very well 
known and hence it has to be modelled, which will introduce some systematic 
uncertainty. However, integrating over the window in hadronic invariant 
masses relevant to b — > u transitions, we have shown that the dependence 
on the shape function is much smaller than the uncertainties induced by the 
quark mass and by the truncation of the perturbative series. 

The the ratio between the semileptonic rates including a cut on the hadronic 
invariant mass and the semileptonic rate without a cut yields U P to 

a quantity r, which we have computed in leading twist approximation and to 
order a s . Based on our calculations the uncertainty in this quantity is 20% 
leaving us with a 10% theoretical uncertainty in the determination of V u b- 
This method is thus one of the cleanest possible to obtain V u b at the ongoing 
5-factory experiments. 
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r = 1.57 ±0.3 



(24) 




(25) 
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Appendix 

We present here the convoluted spectra in some more detail, showing in 
particular a convenient way of calculating the contribution from the terms 
proportional to a s . The convolution of the partonic spectra results in the 
leading-twist approximated rate of B — > X U W decay. The resulting distribu- 
tion can be written in the following form: 

dr r 1 r 2 - x /-zma* i-m b ^YP arton 

dx I dx p / dz / dm*f(m*) — — — x 



dM x dEi Jo Ji-x Jz min J-oo dxdxpdz 

5(M 2 X - m* 2 z - x p A*m* - A* 2 )5(^ - ^). (A.l) 

The partonic rate can be split into the tree level term and the 0(a s ) correc- 
tion, 

^£parton ^^partonfi _|_ ^ppqrfon,! ^ 

3ir 

and the convoluted rate can then be divided up accordingly. Using the delta 
functions, the integrations involved in these rates can be simplified, yielding 



m b /2 f M* ut dr(°)(S -> X u ev e ) 

dEi / dMl \. — - J = 

E cut l Jo x dM x dE t 

rMg 12 rl rx 9 / 9 K, \ f ir partori '° 

/ d El dx dyiff^)*——, (A.3) 

JE cut J2E 1 /m b Jy min x V x J dxdy 

y mm -maxjO, {2Ei/x){Mb _ 2Ei/x) ] (A.4) 
for the tree-level contribution, while the a s correction contributes 

JE cut l Jo X dM x dE t 

rM B /2 A r2-x rz ma x 2 /2i?A dF^* ™' 1 



IE, 



cut 



dE t dx dx p dz-ff—L) A.5 

J2E l /M B Ji-x Jz min x \ x J dxdx p dz 
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where 

f 0, ior x p <l, (A.6) 



Xp — 1 for x p > 1, 



and 



(A.7) 

With these formulae, it is easy to numerically convolute the Born approxi- 
mated rate. However, it is useful to eliminate one integral from the convolu- 
tion of the a s term. To this end, the integral over z in Eq. (A. 5) has been 
performed analytically, so that we use the prime function F(x, x p , z) defined 

as 

rz ^ ^partem,! 

J Zmin Todxdxpdz' 

where T is defined in (9). Then, 

F= iF 1 (z)-F 1 (0) + F 2 (z)+F 3 , x p <l, 

\ F^z) - F x {z min ) + F 2 (z) - F 2 (z rnin ) x p > 1, { ' ' 

where 

Zmin %p 1) (A. 10) 

and 



Fi 



1, At i i /o\2 /.2n i q/ 1 i „ lc,\2 1,2 



\n——[x 2 p - 4(x - 1 + x p /2Y/t 2 } +S(x-1 + x p /2y/t i 

x(36 - 39x p + 12x 2 p - 3x 3 p /4) + - In ^-^(x - 1 + x p /2)[48 - 78x p 

+45a£ - 2\x\j2 + 3xJ/4 + (x - 1 + x p /2)(-42 - 15x p + 27xJ 

-21/8xJ)] + (Ax* + A 2 ) In ^ + A 3 log(l -t) + A 4 , (A.ll) 

F 2 = 12[21n 2 z- (81nx p -7)lnz]t;(x-Xp), (A.12) 
^3 — (^p — {241na;p (—1 + 5v — 4v\nx p ) 

+v[-lQn 2 - 60 - 48Li 2 (l - x p )]} . (A.13) 

In the above formulae, 

Ai = (x p — x){xpZx/2 + 21x p x/2 — 6x p v — 12x p x + 3x p v + 7x p x/4) 
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+(2x p zv - 9x p z/2 - 9/10x p z 2 - x\zv + 7x 2 p z + 12x 2 p v + 15x 2 p /2 
-21x 3 p z/20 + 17x 3 p v/2 - 105/4xJ - 5/4x 4 p v + 53/4x 4 p 
-219/160^), (A.14) 
A 2 = (x p - x)(-3Qx p (l - x p x) + Qx 2 p v + 4x 3 p x + 12x - 96?; In 2 + 48v) 

-2Ax p + 48xJ + 16xju - 60a:J - 2xjv + 26xJ - 12/5xJ, (A. 15) 
A 3 = 2A 2 + 9Q(x p -x)(l + x- x p ) ln(l + i), (A.16) 
A 4 = (x p — x)(5x p zx — 2Azx + 12 zv ) + 32x p 2;t> — Slx p z — 21/hx p z 2 

-4x 2 p zv + 40x 2 p z - 81/20x^ + 36^ + 12z - 6z 2 v + 24z 2 , (A. 17) 

and 

x = l-x, w = 2-x-a;p, t = yjl - 4z/x 2 p . (A.18) 
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